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ABSTRACT 

A  method  for  predicting  the  average  monthly  dispersion 
pattern  of  a  passive  contaminant  point  source  is  developed  utilizing 
current  data  in  near  shore  areas.   The  method  applies  turbulence 
theory  concepts  on  a  monthly  basis  to  obtain  the  dispersion  para- 
meters.  Dispersion  parameters  determined  on  a  monthly  basis  agree 
with  similar  values  determined  on  hourly  observations. 


TWO-DIMENSIONAL  DISPERSION  PLUME  PREDICTIONS 

BASED  ON 

NEAR  SHORE  LAKE  CURRENT  METER  RECORDS 
INTRODUCTION 

Recording  current  meter  data  was  collected  in  the  Nanticoke 
area  on  Lake  Erie  during  1968  as  part  of  a  detailed  water  environment 
study.    Previous  work    applied     standard  time  series  methods  to 
determine  energy  spectra  of  currents,  temperatures  and  water  as  well 
as  the  correlations  of  these  parameters  (Palmer,  1969) »   This  work  thus 
provided  detailed  information  on  the  water  movements  and  temperatures 
in  the  area  as  well  as  the  physical  forces  responsible  for  these.  There 
was  also  a  need  to  develop  the  necessary  dispersion  relationships  for 
passive  contaminants  introduced  into  the  area  to  permit  the  assessment 
of  the  effects  of  a  source  of  pollutants  on  the  local  environment u  This 
paper  outlines  a  method  for  the  prediction  of  dispersion  relationships 
based  upon  recording  meter  data0   The  terms  used  in  the  method  originate 
from  turbulence  studies  and  are  defined  where  necessary.  Intermediate 
steps  have  been  omitted  to  permit  a  concise  pre sentation „ 

DEVELOPMENT  OF  METHOD 

J.  O.  Hinze,  1959,  in  his  book  entitled  "Turbulence",  presented 
various  forms  of  dispersion  equations  for  different  types  of  flow  fields „  All 
equations  contain  terms  representing  turbulence  characteristics,  distances 
from  sources  and  mean  velocities  where  applicable.   The  problem  is  thus 
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to  select  the  appropriate  equations  for  the  Nanticoke  area. 

Measured  velocities  using  current  meters  are  two-dimensional/ 
consequently,  any  relationship  describing  the  dispersion  characteristics 
will  have  to  be  two-dimensional.   This  is  not  a  serious  limitation  since 
the  maximum  depth  in  the  study  area  is  only  13m  and  provided  the  long  term 
distribution  of  a  passive  contaminant  is  considered.   Over  a  long  period  of 
time  the  distribution  with  depth  will  be  small  compared  to  the  horizontal 
(parallel  to  the  water  surface)  distribution .  Energy  spectra  forthe  north-south 
and  east-west  correlations  between  meter  locations  indicated  that  the 
dispersion  characteristics  for  the  components  are  different ,  Therefore, 
the  applicable  expression  must  differentiate  between  the  components. 
One  equation  that  meets  these  requirements  is: 

P£«„x>  f  dt.      S^__  exp  k(C*,- OMji  1 5i 


—  00 

(Hinze,  1959,  pu  327) 

where:       ^         =        distance  from  source  in  north- south  direction 

in  cm  u 

Xj         =        distance  from  source  in  east- west  direction  in 
cm. 

S  =        continuous  point  source  cmVsec« 

variance  of  the  displacement  in  the  east-west 
direction  in  cm^ . 

variance  of  the  displacement  in  the  north- south 
direction  in  cm^ . 

mean  velocity  in  the  east- west  direction  in 
cm/  secu 
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P(Xj,  x2) 


probability  of  finding  a  marked  fluid  particle 
at  a  point      ,  *2  • 


t 


time  in  sees. 


Equation  (1)  is  a  special  solution  of  the  general  equation  (Hinze,  1959, 


This  is  the  classical  diffusion  equation  for  molecules  and  heat  in  four- 
dimensional  space  which  states  that  quantities  move  down  concentration 
gradients. 


function  of  the  space  and  time.   To  obtain  dispersion  plumes  from  the 
two  dimensional  (north- south  and  east-west)  current  meter  records ,  it 
is  necessary  to  assume  that  the  dispersion  is  described  by  the  two 
horizontal  diffusion  coefficients,  and  that  the  coefficients  remain  con- 
stant throughout  the  plume.   The  first  assumption  of  two- dimensionality 
is  not  serious,  provided  that  reasonably  long  time  periods  are  considered 
so  that  the  distribution  with  depth  is  uniform,  or  conversely  that  the 
distances  are  far  enough  from  the  source  such  that  the  depth  scale  of 
13m  is  not  significant.   The  second  assumption  that  the  diffusion  co- 
efficients are  constant  with  distance  from  the  source  can  be  justified 
by  developing  the  coefficients  as  long  time  averages  and  invoking  the 
condition  that  the  velocity  field  is  reasonably  uniform  over  areas  of  two 


p.  312): 


where         is  the  diffusion  coefficient  tensor  of  second  order 0 


so  that  the  diffusion  coefficient  is  a 
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or  three  miles.   The  uniform  velocity  field  assumption  has  been  verified 
at  Nanticoke  by  the  time  series  analysis  of  velocities  which  revealed 
significant  correlations  between  meters  022  and  023  (Palmer,  1969),  It 
was  decided  to  restrict  the  assumption  of  constant  diffusion  coefficient 
to  distances  of  5  x  10^  cm  or  50  integral  scales  (Lumley,  1964). 
For  short  times  (Hinze,  1959,  p. 332) 
where 


*y  =  time 

=      Lagrangian  integral  scale 


(3) 


For  long  times  (Ibid) 


where  <££L  |, 


U*  =    2  £  3C<  .  (4) 

r.  O. 

and      £    =    LL^A_V.   (5) 

\jj  =      fluctuating  velocity  root  mean 
square  in  cm/sec. 
— A~l~      Lagrangian  integral  scale  in  cm. 

Since  the  meters  are  fixed,  the  velocity  data  is  Eulerian  not  Lagrangian 

However,  for  extremely  large  Reynolds  numbers,  it  can  be  shown  by 

dimensional  reasoning  (Lumley,  1964,  p. 34)  that: 

where        =      Eulerian  integral  time  scale 
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The  Reynolds  number  for  Nanticoke  is: 

TRe  -  =  Io^kXO*  (very  large) 

where   {_)    =      mean  velocity  in  cm/sec. 
T>    =     depth  in  cm. 

=      kinematic  viscosity  cm Vsec 0 

By  definition 

^  ptt  u(t')uCt'-t)  .  J^dl  TJtU  

where  =      correlation  coefficients 

=      integral  of  the  correlation  coefficients 
e  with  respect  to  time  in  secsu 

{jl      =      fluctuating  velocity  in  cm/sec. 

The  current  meter  data  was  analyzed  on  an  hourly  basis  to 

minimize  alaising  in  the  energy  spectra 0    Consequently,  correlation 

coefficients  were  developed  for  hourly  time  lags.  The  question  arises 

as  to  whether  hourly  lags  are  appropriate  for  determing  "R,    ?  Hinze 

(1959,  p. 4)  indicates  a  sampling  time  interval: 

T      =      length  scale 
mean  velocity 

Using  the  depth  of  13m: 

T      =      7  minutes 
On  the  other  hand,  Lumley  (1964,  p. 39)  indicates  a  sample  time  interval 
of: 

200  x  integral  scale 
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Okubo  (1967)  estimated  the  integral  scale  to  be  approximately  2m  for 
Lake  Erie: 

T      =      200  x  200    =  220  minutes 
3  x  60 

Since  Hinze' s  intervals  are  basically  laboratory  oriented  while 
Lumley' s  are  meteorologically  oriented  the  differences  are  likely  a  result 
of  the  scale  of  the  processes.   The  Nanticoke  characteristics  would 
obviously  tend  towards  the  meteorological  scale.    Consequently,  it  seems 
reasonable  to  utilize  hourly  lag  intervals  for  the  correlation  coefficients 

To  evaluate  the  integral  in  equation  (6),  the  limits  must  be 
reduced  from  infinity  due  to  noise  and  very  large  scale  motions »  An 
upper  limit  of  50  hours  was  selected  (see  Figsu  2  and  3)  as  covering  the 
time  interval  of  major  interest.,   As  it  is  necessary  to  consider  long  time 
periods  for  the  two  dimensional  assumption,  it  was  decided  to  use  a 
month  as  a  convenient  time  as  the  basis  for  determining  dispersion  patterns. 

Velocities  in  the  Nanticoke  area  were  small  (monthly  averages 
for  Meters  023  and  022  were  1.7  and  4Q7  cm/sec.  respectively)  with 
complete  reversals  of  direction  frequently  occurring  within  2  or  3  days 
(see  Figs0  3a  and  b,  Palmer,  1969)  0   This  is  typical  of  near  shore  areas 
on  lakes  Erie  and  Ontario.    Consequently,  it  would  be  extremely  difficult 
to  consider  dispersion  for  short  time  periods „   A  fundamental  time  period 
of  a  month  permits  the  averaging  of  many  reversal su   Obviously,  if  the 
currents  during  a  month  were  very  persistent  in  one  direction  (with  a 
persistence  factor  of  0o8  or  greater),  it  would  be  possible  to  consider 
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shorter  time  periods  with  accuracy. 

The  monthly  resultant  current  vector  was  used  for  the  mean 
velocity   {J%  and  the  standard  deviations  of  the  monthly  currents  was  used 
for  estimates  of  the  fluctuating  velocities  AX  0   The  corresponding  persis- 
tence factors  provide  an  indication  of  the  probability  of  the  dispersion 
pattern  occurring  in  some  other  direction  during  any  day  of  the  month. 
The  values  thus  derived  are  tabulated  below. 


September  1968 


Meter 

*3 

North 
hours 

"3 

East 
hours 

TV 

North 
cm 

A 

East 
cm 

/ 

North 
cm/  sec. 

LA 

East 
cm/  sec. 

022* 

9.6 

19,5 

8  06  x  103 

1.46  x  104 

1.49 

4.94 

023* 

7.25 

9.5 

3.62  x  103 

4.02  x  104 

0.65 

3.08 

u 

North 

cm/  sec. 

u 

East 
cm/  sec  o 

Persistence 
Factor 

North 
cm2/sec . 

e 

East 
cmV  sec . 

022 

0.25 

2.08 

0o44 

1.28  x  104 

7.19  x  105 

023 

0.14 

lol9 

0.70 

2.36  x  103 

1.25  x  105 

for  meter  locations  see  Fig.  1 


Results  of  Other  Studies 


cm2/  sec. 

Csanady  (1964)  4  x  102  (one  at  2  x  103) 

Okubo  (1967)  3  to  6  x  104 

Noble   (1961)  2U44  x  102 

Palmer   (1968)  0„65  to  1  u3  x  102 
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The  above  results  compare  favourably  with  Okubo,  (1967) 
who  used  drogues  at  similar  distances  offshore  but  are  significantly 
larger  than  the  other  studies,,    However,,  as  the  other  studies  are  either 
surface  dye  dispersion  determinations  or  evaluations  in  deeper  water 
it  is  felt  that  a  comparison  is  not  valid 0   A  very  important  feature  of  the 
Nanticoke  results  is  that  they  portray  the  nearshore  processes  more 
precisely o 

The  dispersion  is  more  intense  nearer  the  shore  (Meter  022) 
and  parallel  (east)  to  the  shore „    Both  the  fluctuating  velocities  and  the 
integral  scales  by  the  preceding  method  are  much  larger  than  those 
determined  by  Okubo  (1967)  who  found  the  fluctuating  velocities  to  be 
0o3  cm/sec,  and  the  integral  scale  one  to  10  metreso   It  is  certainly 
unlikely  that  the  integral  scale  would  be  as  small  as  Okubo  estimates u 
Significant  correlations  exist  between  velocities  at  Meters  022  and  023 
(Palmer,,  1969) »   The  discrepancy  in  the  fluctuating  velocities  is  due  to 
the  different  time  periods  considered,, 

The  main  area  of  concern  for  the  dispersion  plume  pattern  is 
from  4  x  10^  cm  to  5  x  10^  cm  (1/4  miles  to  3  miles)  from  the  source 0 
Considering  that  the  short  time  formula  (equation  (3)  )  is  limited  to  dis- 
tances up  to  approximately  10^  cm  while  the  long  time  formula  (equation 
(4)  )  is  limited  to  distances  beyond  10    cm,  the  area  of  interest  is  some- 
where between  these  limits,, 
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The  short  and  long  time  forms  of  equation  (1)  (Hinze,  1959, 
p. 3 27  and  p. 3 29)  for      =  0  are  plotted  in  Figure  4.   It  can  be  seen 
that  the  area  of  interest  is  centred  about  the  intersection  of  the  short 
and  long  time  forms  of  the  equation.   Thus,  it  is  reasonable  to  use  an 
average  of  the  short  and  long  time  equations  for  this  area.    Short  and 
long  time  equations  for   y    in  two-dimensions  (equations  (3)  and  (4)  ) 
were  computed  for  each  computation  point  on  a  square  grid  (see  Fig.  1) 
and  averaged.   The  values  thus  determined  were  then  substituted  into 
equation  (1)  which  was  then  integrated  numerically  on  the  computer 
until  the  last  step  added  less  than  half  a  per  cent  of  the  total.  (The 
integral  equation  (1)  converges  rapidly  with  time  intervals  of  15  hours) „ 

SAMPLE  CALCULATION 

The  numercial  solution  of  equation  (1)  generates  values  of 

P/S  for  various  points  on  a  two-dimensional  grid  (Fig.  1).   For  a  P  not 

greater  than  1.0  at       =  4  x  10^  cm  and  x.^  ~  0,  the  maximum  point 

5  3 

source  strength  is  6.75  x  10    cm  /sec0  (approximately  24  cfs)  for 
Meter  022  and  "S"  of  2.92  x  105  cm3/sec.  (approximately  10  cfs)  for 
Meter  023.    This  means  that  for  distances  less  than  4  x  10^  cm  from 
the  source  it  must  be  assumed  that  a  uniform  concentration  equivalent 
to  source  strength  exists  since  the  equations  do  not  apply  in  this 
range.   For  larger  flows,  the  area  of  uniform  concentration  would  be 
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greater.   This  is  not  unexpected  since  the  equations  reflect  the  fact 
that  larger  sources  will  change  the  local  mean  velocities   LJ,    in  the 
immediate  area  to  accommodate  the  larger  flows.   The  mean  velocities 
determined  by  current  meter  studies  can  therefore  not  be  used  when  no 
actual  source  is  operative  in  this  area.   The  September  dispersion  plumes 

r  O  5  3 

for  a  source  strength  S  =  6.75  x  10D  cm  /sec.  and  S  =  2o92  x  10  cm/sec. 
at  Meters  022  and  023  respectively  are  plotted  in  Figure  5  as  percentages 
of  the  source  strength.   This  means  that  if  the  source  strength  is  10  ppm, 
the  concentration  on  the  5  0  per  cent  contour  would  be  5  ppm. 

DISCUSSION 

The  preceding  method  for  predicting  the  long  term  dispersion 
plumes  provides  a  comprehensive  portrayal  of  the  nearshore  processes. 
It  permits  differentation  of  the  dispersion  characteristics  in  two- 
dimensions  for  various  locations  offshore.    By  using  the  continuous 
records  from  current  meters,  the  method  considers  the  time  dependant 
and  periodic  nature  of  the  currents  serially.   The  periodic  nature  of  the 
current  regime  necessitates  the  development  of  average  conditions  over 
reasonable  time  periods 0    Consequently,  the  resulting  plumes  (see  Fig. 
5)  do  not  represent  the  dispersion  at  any  particular  time,  but  are  the 
average  of  conditions  occurring  theoughout  the  month.   The  method 
also  predicts  how  far  a  new  source  will  alter  the  current  patterns  in  the 
area. 
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The  assumptions  of  the  method  dictate  the  extent  of  its 
application.    The  expression  is  restricted  to  areas  where  the  Reynold's 
numbers  are  large  and  where  the  assumption  of  constant  long  term 
effective  diffusion  coefficient  is  valid. 
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FIG.  4 

COMPARISON   OF   LONG   AND  SHORT   TIME    DIFFUSION  PROBABILITIES 
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